Higher angular momentum Kondo liquids 



O 

o 

(N 

> 
O 



o 
o 



o 



I 

o 
o 



X 



Pouyan Ghaemi^ and T. Senthil^'^ 

^ Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
"^Center for Condensed Matter Theory, Indian Institute of Science, Bangalore, India 560012 

(Dated: February 5, 2008) 

Conventional heavy Fermi liquid phases of Kondo lattices involve the formation of a "Kondo 
singlet" between the local moments and the conduction electrons. This Kondo singlet is usually 
taken to be in an internal s-wave angular momentum state. Here we explore the possibility of 
phases where the Kondo singlet has internal angular momentum that is d-wave. Such states are 
readily accessed in a slave boson mean field formulation, and are energetically favorable when the 
Kondo interaction is between a local moment and an electron at a nearest neighbor site. The 
properties of the d-wave Kondo liquid are studied. Effective mass and quasiparticle residue show 
large angle dependence on the Fermi surface. Remarkably in certain cases, the quasiparticle residue 
goes to zero at isolated points (in two dimensions) on the Fermi surface. The excitations at these 
points then include a free fractionalized spinon. We also point out the possibility of quantum Hall 
phenomena in two dimensional Kondo insulators, if the Kondo singlet has complex internal angular 
momentum. We suggest that such d-wave Kondo pairing may provide a useful route to thinking 
about correlated Fermi liquids with strong anisotropy along the Fermi surface. 

PACS numbers: 



I. INTRODUCTION 

In many rare earth alloys, a periodic lattice of localized 
spin- 1/2 magnetic moments is coupled through Kondo 
exchange^ to a separate band of conduction electrons^. 
An interesting low temperature metallic phase often de- 
velops at low temperature where the local moments are 
absorbed into the Fermi sea by a lattice analog of the 
Kondo effect. The resulting metallic phase is a Fermi 
liquid albeit with strongly renormalized parameters the 
most dramatic of which is a large quasiparticle mass en- 
hancement (leading to the name 'heavy fermi liquid'). 
Concomittantly the quasiparticle residue at the Fermi 
surface is very small (though non-zero). Perhaps most 
strikingly the Fermi surface is large in the sense that its 
volume satisfies Luttinger's theorem? only if the local 
moments arc included in the count of the electron den- 
sity. This and other universal properties of this Fermi 
liquid are usefully understood in a strong Kondo cou- 
pling picture where each conduction electron is trapped 
into a spin singlet state with a local moment - see Ref. 
i 

The singlet 'molecule' formed out of the local moment 
and conduction electron is usually taken to be in an s- 
wave state with zero internal angular momentum. In this 
paper we explore metallic states where this singlet has 
non-zero internal angular momentum. We show that this 
results in Fermi liquid states that have many unusual and 
interesting properties. For instance such states naturally 
have large anisotropics in the effective mass, quasiparticle 
residue and other properties on moving around the Fermi 
surface. Under certain conditions it is even possible for 
the quasiparticle residue to vanish at isolated points of 
the Fermi surface. The excitations at such points do 
not have electron quantum numbers but can be under- 
stood as a neutral fermionic spin-1/2 'spinon'. Remark- 



ably such a spinon excitation emerges without any asso- 
ciated gauge interaction (unlike the other many familiar 
examples^i^ii) . In the special case where the conduction 
band is half-filled the usual s-wave Kondo singlet forma- 
tion leads to an insulating state (dubbed the Kondo in- 
sulator). If the internal angular momentum is non-zero, 
interesting varieties of Kondo insulators become possi- 
ble. For instance we show that in two dimensions with 
a dx^-y^ + idxy singlet, the Kondo insulator has a non- 
trivial quantized electrical Hall conductivity. 

A partial motivation for studying such higher angu- 
lar momentum singlet formation between local moments 
and conduction electrons comes from observations on 
cuprate high- Tc materials. We should immediately em- 
phasize though that the two band Kondo-lattice mod- 
els are not directly expected to be an appropriate mi- 
croscopic description of the cuprates. Nevertheless it 
seems worthwhile to point out some resemblances be- 
tween the properties of the metallic states described in 
this paper and those seen in the cuprates. Extending 
the ideas of this paper to one-band models appropriate 
to the cuprates may be an interesting direction for fu- 
ture work. For some attempts along similar directions 
see Ref. H. The optimally doped cuprates in the normal 
state are metallic and have a full Fermi surface satisfy- 
ing the usual Luttinger theorem^.. Despite this however 
a true sharp electronic quasiparticle peak does not seem 
to exist. The electron spectrum - measured through pho- 
toemission experiments - shows significant anisotropy on 
moving around the Fermi surface^. In particular the 
energy distribution curve (EDC) is narrower along the 
nodal direction as compared to the antinodal direction. 
This kind of anisotropy in the EDC persists into the un- 
derdoped side when a pscudogap opens along the antin- 
odal direction. It also persists into the superconducting 
state obtained by cooling. Collectively these phenomena 
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have been dubbed the 'nodal- antinodal dichotomy': the 
nodal excitations arc much more quasiparticlc-like than 
the antinodal ones. 

The general message from these results on the cuprates 
for the study of correlated metals is that correlation ef- 
fects may not uniformly affect all parts of the Fermi sur- 
face in a metal. Some portions of the Fermi surface may 
be more susceptible to correlation effects than others. 
This will lead to strong correlation induced anisotropy 
in physical properties around the Fermi surface. The- 
oretically study of such effects is difficult because elec- 
tron correlations are typically most easily handled in real 
space which does not readily distinguish between differ- 
ent parts of the Fermi surface. Interesting recent numer- 
ical calculationsli based on cluster extensions of dynam- 
ical mean field theory on single band Hubbard models in 
two dimensions have found such momentum space differ- 
entiation in the electronic properties. The ideas of the 
present paper may be seen as a step toward incorporat- 
ing momentum space differentiation into simple analytic 
treatments of strong correlation problems. Anisotropic 
spectra at different regions of the Fermi surface occur 
naturally in the metallic states described in the two band 
Kondo lattice models studied in this paper. 

Theoretically higher angular momentum Kondo liquids 
are conveniently accessed through the slave boson mean 
field theoryi^ii^ developed to describe the usual heavy 
fermi liquid state. In this theory the local moment Sr at 
site r is represented in terms of neutral spin- 1/2 fermions 
fra [a =t, i)through 

Sr = \flSfr (1) 

The Kondo singlet formation is then signalled by the de- 
velopment of a non-zero expectation value of the singlet 
hybridization amplitude (cj/r) between the /-fermions 
and conduction electrons Cr- This immediately suggests 
that Kondo singlets of higher angular momentum can be 
described by hybridization amplitudes between at site 
r and c^' at a different site r' . Specifically 

brr' = {cUr') (2) 

may be viewed as the wave function of the Kondo sin- 
glet. Thus by choosing the internal angular momentum 
associated with rotations of the relative coordinate f' — r 
appropriately Kondo liquids with higher angular momen- 
tum may be constructed. In this paper we will exploit 
this strategy to construct mean field descriptions of var- 
ious such Kondo liquid states. Within the mean field de- 
scription, higher angular momentum Kondo liquids cor- 
respond to a particular form of momentum dependence of 
the Kondo hybridization amplitude. Momentum depen- 
dent hybridization amplitudes have been previously con- 
sidered in Ref. flil . Recent optical transport experiments 
on the 1-1-5 materials have also been interpreted in terms 
of momentum dependent hybridization amplitude o^^^^ . 

We begin in Section |TT] with a brief review of the 
usual slave boson mean field theory for a heavy fermi 



liquictii. This mean field theory describes an on-site "s- 
wave" Kondo singlet. We then consider in Section Hill a 
Kondo lattice Hamiltonian where each local moment in- 
teracts through Kondo exchange with a conduction elec- 
tron at a neighboring site. Within mean field theory this 
Hamiltonian is shown to stabilize a d-wave Kondo liquid. 
The properties of this state are then studied in Section 
IIVI Next we modify the Kondo lattice model of Section 
mil by introducing explicit Heisenberg exchange between 
nearest neighbor moments. A mean field treatment of 
this model (Section |V| leads to a metallic state where 
the quasiparticle residue vanishes at isolated points of 
the Fermi surface. As mentioned above this remarkable 
state has spinon excitations at these isolated points with- 
out any associated gauge interactions. Next in Sect ion IVTl 
we consider two dimensional Kondo insulators where the 
Kondo singlet has complex internal angular momentum 
of the form dj;2_y2 +idxy We calculate the Hall conduc- 
tivity and show that it has a non-trivial quantized value. 
Thus this provides an interesting example of a quantum 
Hall state in a Kondo lattice. Section IVIIII has some 
conclusions, various details are in the Appendices. 

II. REVIEW OF MEAN-FIELD THEORY FOR 
KONDO LIQUID: s-WAVE SINGLET 

Consider the Kondo lattice model describing a lattice 
of spin-1/2 local moments St and conduction electrons 
Cia (a =T, i) coupled through Kondo exchange^. 

^ = -E^4c, + ^^S..clac. (3) 

The total number of conduction electrons per site is taken 
to be some fixed value n^. We also specialize to two di- 
mensional systems though extension to three dimension 
is straightforward. A useful mean field treatment is ob- 
tained using a fermionic representation of the local mo- 
ment spin: 

= ]\\h (4) 

together with the constraint ///^ = 1 at each site. After 
some algebra the interaction term ( Jk) reduces to: 

fjciclfi + constant 

i 

In the corresponding partition function, this takes the 
form 

e/rf^E./.c.e./. (5) 

Using a Hubbard-Stratonovich transformation, this can 
be decoupled as 

1 \DV\ exp + ^yi}.r)U^ + r)c-J,) 

(6)' 
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Exact integration over V{i,T) will produce the original 
interaction term. But in mean-field approximation, we 
do the integral using saddle-point method. In addition 
the constraint fjfi = 1 will be implemented on average. 
We will choose the mean-field solution with V{i,T) = V 
which leads to the following mean-field Hamiltonian and 
self consistency equation for V: 



H 



"^tclcj+h.c 

■^f^fflfi + vY,ifhi+4fi) (7) 



V 

iflfi) 



Jk 

2N 
1 



Y.^flci + h.c.) 



(8) 
(9) 



We have introduced a chemical potential term for the 
/-fermions which serves to set their average number per 
site to be one. In addition the Fermi energy Ep of the 
hybridized quasiparticle must be determined by requiring 
that there are a total of 1 -|- ric fermions per site. Only 
states with energy less than Ep are filled in the ground 
state. This mean-field Hamiltonian is quadratic and can 
be diagonalized. To do so, first we write the Hamiltonian 
in the Fourier space: 

H = Y, ek4ck + t^fflf^ + ^ E(/k^k + 4/k) (10) 

fe k k 

with Ek = {cos{kx) + cos(fcj,)). Now it is easy to 
derive quasi-particle dispersion relation for two bands: 



± \ ( 



Ek - flf 



)2 + y2 



(11) 



With this in hand, we can solve the self consistency con- 
ditions and calculate physical properties. For instance 
the density of states at the Fermi energy is exponentially 
large in 1/Jk leading to the large mass enhancement at 
small Jk- 



III. d-WAVE KONDO LIQUID 

In the last section we considered the case where the 
singlets between the local moments and conduction elec- 
trons are formed with zero internal angular momentum 
(i.e. (fjci + h.c.) is constant). When we assume that the 
singlet is formed between the local moments and con- 
duction electron at the same site (i.e. singlet is local), 
the s-wave is the only possibility for internal state of the 
singlet. In this section we consider modifying the Kondo 
lattice Hamiltonian so that Kondo singlets with non-zero 
angular momentum are favored. To that end we consider 
a generalized Kondo Hamiltonian: 



H ■ 



iiJ) 



(12) 



where the Kondo exchange term is not limited to the local 
moments and conduction electrons at the same site. The 
simplest case is for to be non-zero only when i and j 
denote the nearest neighbor sites: 

Jij = Jk (i,j) nearest neighbor (13) 
Jij = otherwise (14) 

With this choice and after some algebra, similar to what 
we did in the last section, the Hamiltonian reduces to the 
following form: 



{ij) (ij) 



(15) 



Now proceed as before by decoupling the interacting part 
of the action using an auxiliary field, which this time lives 
on the bonds of the lattice, instead of the sites: 



/ 



[DV] 



(16) 



As before, to get the mean-field Hamiltonian, we'll do the 
integration over V using saddle point approximation (in 
addition to imposing (f-fi) = 1 on average). There are 
different saddle point solutions (they are basically differ- 
ent local minima in functional space). We will consider 
two sohitions named as s-wave and d-wave (the reason 
for choosing these names will become clear later). In 
s-wave solution we consider V{i,j,T) = V on all bonds 
and in d-wavc, we consider V{i,j, r) = ^ on a; bonds and 
y(i,j, t) = —V on y bonds: 

[ d ) ' 

(17) 

It is obvious from this form that these correspond to 
s-wave and d-wave internal state for the singlet. In mo- 
mentum space s-wave and d-wave Hamiltonian have the 
following forms: 

ff«(k)=^6kctck + 5^M//^/k 
k k 

-H F^(cos(fc^) + cos{ky)){flck + 4/k) 

(18) 

Hd{k) = X]ek4ck + E M//k/k 

fe k 

+ F^(cOs(fc,) - COs(fc,))(/Tck + 4/k) 

(19) 
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With this quadratic mean-field Hamiltonian, we can also Va 

derive the spectrum and the mean-field parameter V self- 

consistently: 




Ei{k) = 



2 



^ekM/)2 + y2(cos(fc,) - cos(fc,))2 



(21) 



The ground state is obtained by filling all states upto 
the Fermi level. As before the Fermi energy is fixed by 
requiring that there are I + Uc fermions in the ground 
state per site. Note that in both s and d wave cases the 
-|- band lies entirely above the — band (i.e min{E^) > 
max{E_)). Therefore in the ground state only the E_ 
levels are occupied. The ground state energy is 



E- 



gd 



s!_''*(k))£;!_^''(k) 



(22) 



The self-consistency equations that determine V and /z/ 
are now readily obtained. For instance we have 



2N 



Y^eiE 



F - E'j^'ikyjx 

{cos{kx) ± cos(fcj^))^ 

^i^)2 + l/2(cos(fc^) ± C0s(fcj,))2 



(23) 



In equation!^ +(^) corresponds to s- wave (d- wave). In 
addition we need to impose the condition (///i) — 1 and 



Let us first specialize to the half-filled case ric = I- 
In this case the microscopic model has a particle-hole 
symmetry under which 



y t 



(24) 



where = (— l)(^+i') = ±1 on the A and B sublattices of 
the two dimensional square lattice. As the total number 
of fermions per site l + ric = 2 in this case, all the iJ_ lev- 
els are filled. At the level of the approximate mean field 
Hamiltonians Eqn. 1171 under the particle-hole transfor- 
mation 



fia 

V 



V 



(25) 
(26) 
(27) 



Thus a particle-hole symmetric mean field state (which 
we assume) requires /i/ = 0. 



FIG. 1: V versus J for s-wave and d-wave self consistency 
equation 



With this in hand we proceed to compare the ground 
state energy of the s-wave and d-wave mean field Hamil- 
tonians in half-filled case. First consider self-consistency 
equation for s-wave case: 



1 



1 

2N 



E 



{cos{kx) + cos{ky))^ 



v/(f )' + ^'(cos(A:,) + cos(fc,))2 



(28) 



It is obvious that the right hand side of equation [55] is 
monotonically decreasing function of \V\. So it has its 
maximum value for V — 0: 



1 



Jk 2NV2i 



I cos(fc^) -t- cos(fcj,)| 



(29) 



The right hand-site of the above equation is finite. So 
there is a maximum value of 1/Jfc for which we can find 
a solution for V in equation [28l On the other hand the 
self consistency equation for d-wave at /x/ =0 takes the 
form: 



1 



1 

2N 



{cos{kx) — cos(fcj,))^ 



^{fr + V^cosik^) - cosiky))-^ 



(30) 



Like s-wave, the right hand site of equation is maxi- 
mum for y = 0. But in this case the right hand side 
of equation [3D] is divergent as ^ 0. So for d-wave, 
the self-consistency solution exist for any value of J^. A 
schematic graph of V versus Jk for s-wave and d-wave is 
plotted in figure [1] 

As it is also clear from the figure, there is at least a 
region of small Jk for which V for d-wave is larger than 
V for s-wave (in fact we can have a region with finite V 
for d-wave but s-wave V as small as we want). Now if we 
look at the lower band spectrum (eqn. [^D] and [5T|l this 
proves that there is at least a region for small Jk where 
d-wave ground state energy is lower than s-wave. 

Upon moving away from half-filling, we expect that 
the value of ground state energy starts changing contin- 
uously; so we still have at least a region of small Jk and 
small doping, for which d-wave mean-field Hamiltonian is 
a better approximation than s-wave. We have confirmed 
this by a direct numerical solution of the self-consistency 
equations. 
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IV. PROPERTIES OF d-WAVE KONDO LIQUID These results were also checked and confirmed numeri- 



cally. 



We now study the properties of the d-wave mean-field 
state away from half-filling. We will see that d-wave 
singlet formation provides a natural route toward very 
anisotropic properties over the Fermi surface. Before 
getting into the details of these properties, we will ex- 
amine some aspects of d-wave dispersion in equation 1211 
A schematic graph of Fermi surface for finite doping in 
the first Brillouin zone, as well as dispersion along a spe- 
cific path, are given in figures [2] and [3] respectively. One 
important feature of the spectrum in figure [3] (which is 
also easily proved analytically in appendix [X| is that the 
maximum value of energy in the Brillouin zone is fj,f, 
that is the energy of the points along diagonal direction, 
where fij < e^. This shows that Fermi surface crosses 
the diagonal direction at a point where < (since 
Ef has to be less than /i/). In diagonal direction the 
spectrum is E^ = '^^^^^^ — \^hzJiL\^ which can be written 
as: 

^^k = for ek < /!/ 
= A*/ for ek > M/ 

combining this with the condition E^ < fJ-f shows that 




FIG. 2; Fermi surface in the first Brillouin zone. Occupied 
region is plotted in gray. 

Fermi surface should cross the diagonal direction at a 
point where the quasiparticles energy is equal to ek- This 
result is in fact the basic foundation for the interest- 
ing properties that will be studied in the next sections 
i.e. over the Fermi surface, in diagonal direction, quasi- 
particles are c-electron type; the f-electron properties ap- 
pear, as we go away from the diagonal direction. In the 
following section we will present the result of numerical 
calculation of properties over Fermi surface. 

The self consistency equation for parameter V could 
also be studied analytically (appendix IB|) . We see that 

for d-wave Kondo model, similar to the on site s-wave 

c 

Kondo, we have V oc e -'k . Density of states at Fermi 
energy is also studied analytically (appendix [C|, it again 
shows exponential dependence on coupling constant: 

p{Ef) oc V cx e 




FIG. 3: Spectrum along the lines (0, 0) (n, n) (0, vr) 
(0,0) 

In the following sections we present the effective mass 
and quasi-particle residue on fermi surface, as a function 
of angle (from the center of Brillouin zone). We will see 
that these properties are very anisotropic as we move 
along the Fermi surface. 



A. Effective mass 

The effective mass of quasi-particles is defined through 
second derivative of energy with respect to k± which is 
the momentum in direction perpendicular to Fermi sur- 
face: 



1 

m 



d'^Ej 



(32) 



A closed form for the effective mass could be derived. In 
figure m we have plotted the inverse effective mass over 
Fermi surface near diagonal direction. As we expected, 
along diagonal direction the effective mass matches the 
electron effective mass. As we go far away from the diag- 
onal direction, effective mass becomes very large. This is 
because away from the diagonal direction, the quasiparti- 
cle is essentially an /-fermion with some weak admixture 
with the c-electron. Between these two limits we see the 
strange anisotropic behavior where second derivative of 
energy goes from electron type, positive value to large 
negative value. A comparison with quasi-particle residue 
plot (figure ini) shows that this behavior occurs where the 
quasi-particles are a complete mixture of c-electron and 
/-fermion. 

This weird behavior could be traced by looking closer 
at the spectrum in the Brillouin zone. First consider di- 
agonal direction (figure[5]). We see that moving along this 
line we go from region where E^ = Ck to region where 
E^ = M/ (see figure [3]). So if we plot first derivative of 
(31) energy with respect to k±_ (which is diagonal direction 
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44.2 44.6 45 45.4 
Angle in BZ 



45.8 



FIG. 4: Second derivative of energy respect to k± in direction 
perpendicular to fermi surface. 



for points along diagonal direction), as seen in figure O 
there is a jump from a finite value to zero at the point 
where ek = /i/. In this direction, as shown before, at 
the Fermi surface = ek so there is a finite value for 
second derivative. As we move away from diagonal direc- 
tion, the jump softens slightly (because Vk moves slightly 
from zero) . But still first derivative have a large decrease 
in a small interval and so the second derivative is large 
negative number; interestingly, Fermi surface does cross 
this region at some points. This behavior is very strange. 
In fact we see points at which the effective mass is much 
smaller than free electron mass. 



45,00 degree 44,86 degree 44.72 degree 



dE/dk 



FIG. 5; First derivative of the energy with respect to k±. Ver- 
tical axis, in each plot, is g^'j^ and k± is the momentum in 
direction perpendicular to the Fermi surface at the point de- 
fined by the angle on the top. Vertical line shows the position 
of the Fermi surface 



B. Quasi-particle residue 



of which the Hamiltonian is diagonal. The c-electron 
annihilation operators could be written as: 



Ck 



Uk 7k 



Wk 7k 



(34) 



Mk 



VkWk 




44.2 44.6 45 45.4 45.8 
Angle In BZ 

FIG. 6: Quasi-particle residue on the fermi surface 

Using this form (and the fact that Hamiltonian is di- 
agonal in terms of 7^ operators) : 

(c(k,T)ct(k,0)) = u^(7,^(r)7,[^(0)) +^;^(7k (^)7k^(0)) 

Rewriting the Hamiltonian in term of 7 operators {H = 
Ek(7k Sk + 7k Sk)) we get: 

(7^(r)7k^(0))=e-^^ 
Putting these into equation [551 gives: 



G{k, iu^) 



lUJu 



(35) 



By analytically continuing this to real frequencies and 
taking the imaginary part we get the spectral functioni^: 

A{k,uj)^ul6{iu-E+)+vld{uj-E^) (36) 



The electron-electron green function is given byi^: 
G{k,iuj^)^ dTe"^"^Trc{k,T)c^{k,0)) (33) 



Suppose 7j^ and 7j^ are annihilation operators for quasi- 
particles in upper and lower bands respectively, in terms 



This consists of two peaks and the weight under the 
low energy peak will give us the quasiparticle residue 
(Z). We have plotted Z for the points near diagonal di- 
rection in figure [6l Again as we expect, the quasiparticle 
residue near diagonal direction is of order one (electron 
type excitations) and gets very small away from diagonal 
direction. 
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V. d-WAVE KONDO LIQUID IN A 
KONDO-HEISENBERG MODEL: FERMI LIQUID 
WITH SPINONS 

We now study the properties of the d-wave Kondo liq- 
uid state in a model which allows for explicit Heisen- 
berg exchange interactions between the local moments. 
Remarkably we will show that the quasiparticle residue 
vanishes at isolated points of the Fermi surface in such 
a state. The excitation at such points is a free neutral 
fermionic spinon. In a subsequent Section, we show that 
this spinon survives even when fluctuations beyond the 
mean field are included. Thus this d-wave Kondo liquid 
is a Fermi liquid state that supports spinous at isolated 
Fermi points. Specifically we consider the model 

H = '^tijc\cj + JK^Si.Sj + JH^Si.Sj (37) 

(ij) (ij) (ij) 

We proceed as before using the d-wave mean-field ap- 
proximation to treat Si.Sj term; but here, we have an- 
other interacting term, which is the direct Heisenberg 
exchange between the local moments. Expressing this 
in terms of the /-fermions gives rise to a four fermion 
term. We treat this in mean field theory as well. While a 
number of different mean field decouplings are possible, 
we focus here on one in the particle-hole channel which 
endows the /-fermions with a uniform non-zero hopping 
X- This is in turn determined self-consistently through 
the equation 

X^Jniflf,) (38) 
Using this, we get the following mean-field Hamiltonian: 

k k 

+ V^{cos{kx) - cos(fcj^))(/i[ck -f 4/k) 

(39) 

where e/k = A*/ — X {cos{kx) + cos(fcy)). With this 
quadratic Hamiltonian, we can again get the spectrum: 



identified with one of the bands). Consider the quasi- 
particle residue on each band. From Eqn. 1361 it is clear 
that on the Fermi surface of the band, quasi-particle 
residue is given by while on the other sheet (associ- 
ated with E^) it is given by u^. and Uk are defined 
in Eqn. [M] and satisfy 

V{cos{kx) — cos(fcy)) 



ul + vl = l (41) 

Using this we see that the Fermi surface of the — band 
has large quasiparticle residue and thus has essentially c- 
electron character (with weak admixture to /-fermions). 
On the other hand the -I- Fermi surface has small quasi- 
particle residue and has essentially /-fermion character 
with weak admixture to c-electrons. Following reference 
(IS] we name the — and -I- Fermi surfaces as cold and 
hot surface respectively. Remarkably the quasiparticle 
residue on the hot Fermi surface vanishes at four isolated 
points (which are along the diagonal directions). At these 
four points the excitation is a pure /-fermion with no ad- 
mixture to the c-electron. Thus at these isolated points 
the excitation is a neutral fermionic spinon even though 
spinous do not exist elsewhere on the Fermi surface. 




Cold Fermi surface 

FIG. 7: Fermi surface of the d-wave Kondo liquid in the 
Kondo-Heisenberg model 



i^f (k) ^'J^± 



(40) 



Now let us have a closer look at this Hamiltonian and 
spectrum. When U = 0, local moments and conduction 
electrons are not hybridized and we have two separated 
Fermi surfaces. The electron Fermi surface is identified 
by spectrum ek and is small. Spinon Fermi surface, con- 
tains one moment per site and covers half the Brillouin 



zon( 



.19 



Now assume turning on non-zero V . If V is small 
enough, both bands intersect the Fermi energy. The re- 
sulting Fermi surface then consists of two sheets (each 



VI. QUANTUM HALL KONDO INSULATORS 

In this Section we describe an interesting Kondo insu- 
lating state that is possible if the Kondo singlet has non- 
trivial internal angular momentum. Consider the gener- 
alized Kondo Hamiltonian again p2p and assume Jij is 
non-zero for nearest ( Ji) and next nearest neighbors (0/2)- 
Now, as before, we decouple the interacting term using 
an auxiliary field and consider the saddle point solution 
with nearest neighbor V{T,i,j) the same as d-wave, and 
next nearest neighbor equal to iV2 along one diagonal 
direction and —iV2 along the other. Clearly this core- 
sponds to a Kondo singlet with internal angular momen- 



tum dr^ 



idxy , and describes a mean field state that 
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spontaneously breaks time reversal symmetry. This leads 
to the following mean-field Hamiltonian: 

k k 

+ 'Yicos{k^) ~ cos{ky)){flck + 4/k) 

k 

+ ^2 ^ j sin(fc^) sin(fcy)(4/k - /^Ck) 

(42) 

We now specialize to a half-filled conduction band Uc — I. 
Diagonalizing this Hamiltonian, we readily see that the 
ground state in this case is a Kondo insulator with a gap. 
However we now show that the broken time reversal sym- 
metry leads to a non-vanishing quantized electrical Hall 
conductivity. Thus this state provides an interesting ex- 
ample of a quantum Hall effect in a local moment system 
driven by the Kondo effect. To expose this physics it 
is convenient to rewrite the Hamiltonian in terms of a 
two-component fermion operator i/ik^ 



V'k 



Ck 

/k 



We have 



Ek + 



+ m(k).T I ip[ 



(43) 



Here r are Pauli matrices operating on the two compo- 
nents of V'k and mk is a 3-dimensional vector defined in 
the two dimensional Brillouin zone: 



mk = Vi(cos(fca;) — cos{ky))x 
+ V2 sin{k^) sin{ky)y 



(44) 



Note that if |m| ^ over the Brillouin zone, TO/|m| is 
topologically a mapping from a torus (Brillouin zone) to 
the unite sphere. To calculate the Hall conductivity, we 
need the current operator. Note that in this mean field 
Hamiltonian, charge conservation symmetry is realized 
as invariance under i/jk e^^/^k- Thus though the /- 
fermions start off as neutral, the mean field condensation 
of V has endowed them with physical charge. This ob- 
servation then leads to the current operator: 



d 



H 



(45) 



smoothly to a mapping for which integration |46] could 
be calculated easily and gives twice the surface area of a 
unit sphere so that 



2e^ 

h 



(47) 



This is also confirmed with direct numerical integration. 
The essential physics is also simply illustrated by the fol- 
lowing argument. As the relevant integral is a topological 
invariant we first imagine a smooth deformation to make 
V2 infinitesimally small. In the limit when V2 = the 
gap closes and |TO(k)| has four zeroes in the Brillouin 
zone on the diagonal directions where Ck = A*/. At low 
energies, the physics will be dominated by modes near 
these nodes. On turning on a small value of V2, we ex- 
pect that the universal physics is still correctly captured, 
in an approximation that is legitimate for modes near the 
nodes. So we expand m(k) near these nodes and take the 
integral over k^ and ky from —00 to 00. After the expan- 
sion, we get (for the node at (fep, kp) in the first quadrant 
of the Brilloun zone): 



Vi sm{kp) 
+ V2 sin' 
+ tsm{kp) ^ 



X 

\kp)y 

— z 



(48) 



Here 



ki — kp and cos(fcp) 



4t 



After renaming 



as V, 



Qy as kx, Qx + Qy as ky, T^isin(A:p) as v^, tsmkp/2 
and V2 sin(fcp)^ as A, along with 90° rotation of m 



around x axis, we get the following simple form: 



m(k) = VxkxX + Vykyij + Az 



m{k) 



(49) 



points 



From this we trivially see that TO(k) — |^^^..^| 
along z at kx — ky — and points along radial direc- 
tion in kx,ky plane as ^/c^ -I- ky —> 00. So this vector 

covers half the sphere and contributes 27r to the integral 
in Eqn. 1461 Exactly the same contribution arises from 
the other nodes and so we get 4 x 27r = Btt. Putting this 
in l46l gives the result quoted above for the electrical Hall 
conductivity. 

Thus the dx2-y2 -\-idxy Kondo insulator has a quantized 
electrical Hall conductivity. A similar result - a spin and 
thermal quantum Hall effect - was established for two 
dimensional dx^-y^ + idxy superconductors in Ref. [2^ . 



We can now calculate axy using the Kubo formula^°. The 
details of the calculation are in Appendix [Dl We get 

oxy = / d-k'^t^^^,^ (46) 

The value of axy is then invariant under smooth change 
of rfi. It is thus a topological invariant of the kind con- 
sidered previously by Volovik^^. We can now change m 



VII. BEYOND MEAN FIELD THEORY 

Through out this paper we have treated the Kondo 
and other interactions within the slave boson mean field 
theory. Now we consider the role of fluctuations beyond 
this mean field theory. It has been known for a long time 
that the important fluctuations in such mean fleld theo- 
ries are gauge fluctuations^^. Specifically the redundant 
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representation of the local moment spin Si in terms of /- 
fermions leads to invariance of the resulting theory under 
local U{1) gauge transformations^^ where 



The hybridization field in turn transforms as 



(50) 



(51) 



Note that the hybridization field also carries physical 
electric charge equal to the electron charge. As in the 
usual on-site slave boson theory, the mean field state is a 
Higgs phase where the hybridization field is condensed. 
Consequently the gauge fluctuations acquire a gap and 
can be integrated out. The Higgs condensate also im- 
plies that the U{1) gauge charge of a localized /-fermion 
will be screened to produce a gauge neutral object that 
carries physical electric charge 1. Thus exactly as in the 
usual on-site slave boson theory, the gauge fluctuations 
will play an innocuous role in these mean field states, and 
the mean field theory provides a correct qualitative de- 
scription of the universal features of these phases. These 
considerations apply even to the d-wave Kondo liquid of 
Section FVl where the hybridization amplitude vanishes at 
four points of the 'hot' Fermi surface. Thus the vanishing 
quasiparticle residue at those four points (and the asso- 
ciated spinon excitations) survive beyond the mean field 
approximation. 



VIII. DISCUSSION 

The considerations of this paper may be useful in fu- 
ture theoretical work on a number of different correlated 
electron materials. The example of the cuprate materials 
shows that correlation effects may not uniformly affect all 
regions of the Fermi surface. The result is strong corre- 
lation induced anisotropy along the Fermi surface. The- 
oretical approaches to addressing such effects are ham- 
pered by the difficulty that correlation effects are easiest 
to handle in real space and not in momentum space. In 
this paper in the specific case of Kondo lattices, we have 
shown how to incorporate momentum space information 
into the Kondo singlet formation that determines the fate 
of the local moments at low temperature. We explored 
the properties of metallic 'Fermi liquid' states driven by 
Kondo singlet formation in a channel with non-zero in- 
ternal angular momentum. Through out the paper we 
focused on two dimensional systems, though extension of 
our results to three dimensions is straightforward. We 
showed that such metallic states naturally have strong 
anisotropy of the quasiparticle effective mass and residue 
on moving around the Fermi surface. In some cases the 
quasiparticle residue even vanishes at four isolated points 
on the Fermi surface. The excitations at such points may 
be thought of as neutral spin-1/2 spinous that occur with- 
out any residual 'gauge' interactions. Thus these states 
provide interesting examples where strongly anisotropic 



quasiparticle residues are naturally built into the symme- 
try of the states. We also studied two dimensional Kondo 
insulators driven by Kondo singlet formation with com- 
plex internal angular momentum and showed that they 
have a quantized non-trivial electrical Hall conductivity. 
Such Kondo insulators thus present an interesting sit- 
uation where a quantum Hall effect occurs due to the 
Kondo effect. Exploiting the ideas of this paper to de- 
velop techniques for thinking about the angle dependence 
of correlation effects in momentum space is an interesting 
challenge for the future. 
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APPENDIX A: MAXIMUM ENERGY IN 
BRILLOUIN ZONE 



As derived before, lower band energy at each point of 
Brillouin zone (BZ) is given by 



(Al) 



where Vk — V{cos{kx) — cos(fcy)). So Vk = only along 
diagonal directions. From Ell we see for all points in BZ: 



K < 



(A2) 



Right hand side of the above equation is equal to /ij for 
^if < eia and ek for et < M/- Putting these together we 
get: 



r -^k ^ '^k for ek < M/ 
1 < M/ for ek > M/ 



(A3) 



You see that all points in BZ, obviously, satisfy one the 
two conditions (ek < M/ or ek > /i/). So we get for any 
point in the BZ: 



(A4) 



A closer look at equations lAll and IA21 shows that equal- 
ity in equation IA41 could be only for the points along 
diagonal direction with ek > /i/. Such points cover a 
region with zero volume in the BZ. So that for any finite 
doping Ef < fif (see section |IV|) . 
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APPENDIX B: SELF CONSISTENCY EQUATION Using this for the result of integral we get: 



Analytic treatment of self consistency relation is easier 
and much more clear at zero doping {^f = 0). So we start 
with this case and later, we can argue if our result might 
be modified at finite, but small, doping. In continimm 
limit, self consistency relation (at zero doping) have the 
form: 



1 



1 

2n 



(fk 



{cos{kx) — cos{ky)Y 



(27r)2 ^(f)^ + y^(cos(fcJ-cos(A;,))^ 



(Bl) 

where n is density of lattice sites. We saw before that, 
this equation for V has a solution, no matter how small 
Jk is (see section ??). This was because, the integral 
is divergent &t V = Q. The divergence comes from the 
region where Ck is small. So to study the behavior of this 
integral, it is enough to look at the points where Ck is 
small. Because of the symmetry of BZ with respect to 90° 
rotation, it is enough to consider region with kx,ky > 0. 
The point with ek = in this region, are points with 
kx = TT — ky. We want to look around such points so if we 
change the variables to ky ky and k^ ^ q = ky+k^ — n, 
we just want to look at region with small q. Putting these 
variables in to in to lBll and expanding up to first non-zero 
order in q, we get: 



1 



1 

2n 



dkydq 



{sm{ky)q + 2cos{ky))'^ 



,y{tsm{ky)q)^ + V^{sm{ky)q + 2cos{ky))^ 

(B2) 

As discussed before, we like to study the singular behav- 
ior around q = so we can ignore the q dependence in 
non-singular terms. The relation simplifies to: 



1 



1 

2n 



dkydq 



{2cos{ky))'^ 



(2^)' v/(isin(fc^)(7)2 + F2(2cos(fcy))2 



(B3) 

Now we first perform the integral over ky . Note that the 
limits for this is q dependent (i.e. 0+O{q^) and 7r-|-0(q^)) 
but these also could be ignored since they have no effect 
on singular behavior: 

1 _ 1 f dq r dky {2cos{ky))^ 
Tk~2^j2^Jo (27r) ^{ts\n{ky)qY + V^{2cos{ky)f 

Integral over ky could be perform exactly which gives: 



-,F,[-,-,2,l~^) 



(B5) 



here, 2F1 is Hypergeometric Function. We are in inter- 
ested in small q behavior so 1 — « — The expres- 
sion could be simplified using the idcntitj*2i: 

2i^i(a,6;c;z) = (1 - z)"" 2^1 (a, c - 6; c; 



1 I f dq 2 ,11 , , , 

/ -i^_= 2F1 -,-;2;l B6 
Jfe 2n J 27r _|_ y2 2 2 



where we used the fact 



is convergentss we get: 



1. Since 2Fi(i,i;2;l 



Jk 



A 



dq 



^/q^+V^ 



(B7) 



where A is a finite constant. Integral over q is now trivial 
and leads to result mentioned before (see section [TV]) : 



y (X e -'h 



(B8) 



After dopinping, continuum form of self consistency re- 
lation changes: 



1 



1 



Jk 2n J (27r)2 ^ 



{coB{kx) ~ cos{ky)y 



)2 + l/2(cos(/c^) -cos(A:j,))2 
(B9) 



Again dominant contribution comes from the region 
where |ek — M/| is small. Similar expansion could be car- 
ried out but this time the points we look at are close to 
different curve (defined by cos(A;2;) + cos(fcy) — ^ which 
is slightly away from kx — tt — ky for small fj,f). We 
don't expect (and in fact the resulting integral shows) 
that there is not much change in singular behavior as V 
goes to zero. So we expect that behavior seen in lBSI holds 
and in fact this is confirmed with numerical studies. 



APPENDIX C: DENSITY OF STATES 

To get density of states we use the general formulap^^: 

"^^'^^ Jwf^^^'-^^^ ^^^^ 

Again symmetries in BZ are helpful. First of all we 
can do the calculation on one patch of Fermi surface (in 
kx < and ky < quarter). Also because of reflection 
symmetry respect to diagonal, we can do the calculation 
for the points on Fermi surface which are above the di- 
agonal and double the result. Now in this region we can 
safely (since Ck is one to one function of {kx,ky)) do the 
change of variable, and work with ky and Ck: 



p{Ef) 



dei. 



(27r)2 2V1 - (t - cos(fc,)) 



(C2) 
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Note that ^\ — — cos(fcj,))^ — sin(fc2;) is never zero 
in the part of BZ we are integrating over. In terms of 
new variables: 



)2+y2(|_2cOs(fc,))2 

(C3) 

Using this we first do the integration over which leads 
to: 



1 



(27r)2 2^1 - (i - cos(fc^))2 



where m(k) is defined in and a(k) = ^^^^^iil Also 
is the shorthand for Putting this form in Kubo 
formula!^ we get: 



/■ e~*"^ - 1 
J ^ 

^(V'^(r)(a^a(k) + 9^m(k).<7)V(r) (D2) 



/c,/c' 



V'^(0)(9,a(A;') + dym(k').a)i}{W)) 



4 



+ (yV2t)(§-2cos(fcy)) 



^2(S^_2cOs(fc^))2 



(C4) 



here, e is solution of equation: 



^k - 

for Ek, which is a simple second order equation but we 
avoid presenting the answer which is not necessary for the 
rest of calculation. Now to proceed we need to divide 
to points in three different regions. The first region is 
defined by the points where: 



-Ck 



for these points contribution to IC4I is of order one. 
The other regions are where: 



Ek 



2cos(fcy))2 



(C5) 



This contains most of the points on the Fermi surface, 
since from IB8I we know that V is small. Now we expand 
expression in prentices in IC4I up to first non-zero order 
in V'^: 



PiEf) 



T/2 

^(1 - Sign{ei, - /z/)) + -^/(e, h) 



(C6) 

here /(e, ky) is convergent function. Wc don't need the 
detailed form of this function though. The important 
property of this integral is that for points with ek < A*/ 
(which corresponds to points near diagonal as discussed 
in |X| contribution to integral is of order one. We ex- 
pected this since in this region quasi-particles are more 
free electron like. However for points with ek > (which 
are away from diagonal) contribution is large and propor- 
tional to V~'^. Putting all these together we get a large 
density of state: 



2C 

p{Ef) cx e -'k 
APPENDIX D: KUBO CALCULATION 



(C7) 



Using equation 1451 we can write as: 
J,{^) = ^^^(T)(9^a(k) -f 9^mrk).a)V(r) (Dl) 



Using the Hamiltonian given in 33] we get the green func- 
tion for -0 fields: 



G{iuj, k) = 



1 



iuj — a(k) — m(k).(7 



(D3) 



With this in hand, the Kubo equation reduces to the 
following form: 



1, 'J 



tr 



(9a;a(k) -I- (92,m(k).(T)(9ya(k) -I- 9yTO(k).(T) 
{iE-a{k) - m{k).a) 
1 1 



{i{LU + E)- a(k) - TO(k).(7) {iE - a(k) - ■m{\<i).a) 

(D4) 



We are interested in the limit of above equation as ^ 
0. In this limit we get the following expression for Hall 
conductivity: 



(Pk dE 
(27r)2 2^ 



tr 



{dxa{\<i) -f dxm(k).a) 



(5ya(k) -I- i9j^m(k).(7) 



(iS-a(k) -TO(k).a)2 
1 



{iE - a(k) - m(k).CT) 



(D5) 



Expanding out the sums, we have several terms but tak- 
ing the trace, cancel some of the terms. After integration 
over the E and dropping the terms which are zero under 
the trace gives: 

_ . f d^k 1 
''-y--'J (27r)2 8|™(k)|3 

tr {dxrn{k) [m(k).(T, 9ya(k) + dyrn{k) .a]) 

(D6) 

After doing some algebra on the Pauli matrices and tak- 
ing the trace, we get the relation given in 1461 
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